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Abstract 

We derive the power spectrum V(k) of the density perturbations produced during 
inflation up to second-order corrections in the standard slow-roll approximation for an 
inflaton with more than one degree of freedom. We also present the spectral index n up 
to first-order corrections including previously missing terms, and the running dn/dhx k 
to leading order. 
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1 Introduction 



Inflation generates the primordial perturbations which are the origin of the rich structures, 
such as galaxies and clusters of galaxy clusters, observed in the universe today. The power 
spectrum of these perturbations is constrained to be approximately scale invariant by a 
combination of the measurements of anisotropies on large and small angular scales, as well 
as galaxy surveys M. Moreover, forthcoming measurements such as the MAP and Planck 
satellites and the Sloan Digital Sky Survey will probe the power spectrum with greater 
accuracy, lowering the observational errors further. Planck, for example, is expected to be 
able to measure the power spectrum of the cosmic microwave background with an error of 
a few percent ||. Thus, it is very important to calculate the power spectrum precisely, so 
that we can make use of the observations fully. 

The calculation of the density perturbations for an inflaton with a single degree of freedom 
has been studied for a long time 0, [|, [J, and there have been extensive works on the multi- 
component inflaton case, especially recently f7|, D, ||, [10| . Among these results, those with 



first-order corrections P, |9j are expected to have errors of the order of one percent or smaller 



for the power spectrum |pl]| , which is comparable with the expected errors in the planned 
observations. Thus we should reduce the errors and improve the accuracy of the calculations. 
However, the only existing result including second-order corrections || considered only the 
single-component case. Therefore, it is needed to consider the multi-component inflaton with 
higher order corrections. 

In this paper, we extend our previous formalism || to the multi-component case, and 
calculate the power spectrum up to second-order corrections in the standard slow-roll ex- 
pansionQ We also present the spectral index n up to first-order corrections with new terms 
inconsistently neglected in the previous result ||, and the running dn/d In k to leading order 
for the first time. It is straightforward to check that our results reduce to those of Ref. || 
in the single component case. 



2 The Spectrum for a Multi-Component Inflaton 

The power spectrum V(k) is defined by 

^P(fc)5( 3 )(k-l) = (^ c (k)Kt(l)), (1) 

where ( ) denotes the vacuum expectation value and TZ C is the intrinsic curvature perturbation 
of the comoving hypersurfaces. In the single-component case, the standard result for the 
power spectrum of the curvature perturbation in the slow-roll approximation is 




The multi-component case is more subtle. Detailed arguments are given in ||; here we 
just give the principal idea. Let N = J Hdt be the number of e-folds of expansion. Then, 

1 For a more general slow-roll expansion, see Ref. p2[. 
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on super-horizon scales, we obtain || 

ATL = 5N. (3) 

That is, the change in the curvature TZ between the initial and final hypersurfaces is given 
by the perturbation in the number of e-folds of expansion. In particular, if we choose the 
initial hypersurface to be flat, i.e. TZ(ti) = 0, and the final one to be comoving, we get 

H c (t 2 ) = 5N. (4) 

We take t\ to be some time during slow-roll inflation, and ti to be some time after inflation 
when TZ C has become constant. The relevant scale is assumed to be still outside the horizon 
at t2- Now, N depends on both 0(ii) and <p a (ti), but as t\ is during slow-roll inflation, it will 
be convenient to express this dependence in terms of 4>(ti) and (fi±(ti), where a = <p a — (f)f T 
is the deviation of a from its value <j) a T on the exact slow-roll trajectory passing through 0. 
Then we can writef] 

ON dN ■ 

SN= w mti)+ wJ 4 ' 1{h) - (5) 

However, 5<pj_ is negligible on super-horizon scales during slow-roll inflation. Therefore we 
get 

n c {t 2 ) = 5N = f^<W0 ■ (6) 

Using Fourier expansion and substituting into Eq. ([!]), we obtain the power spectrum for a 
multi-component inflaton 



k 3 v ' v ' d<p a d<f) h 
Now, to leading order in the slow-roll approximation, 

(50 a (k)50 bt (l))^ ^^(k-l)/^ 

Therefore 

( H\ 2 iah dN dN 
V(k) ~ I — h — — — — j- 



aH=k 



(9) 



aH=k 

Defining to be the unit covector in the direction dN/d<fi a and using 
we can write this as 

V{k) ~ ( — ] I — ! . i l l ; 



V a r 7 aH=k 

and so in the single- component case we recover Eq. (fj), but with some greater insight into 
its meaning. 



2 We use the abstract index notation jl3| . A boldface superscript index denotes a vector in the scalar field 
space, a boldface subscript index denotes a covector, contractions are denoted by repeated indices, and the 
metric /i a b and its inverse /i ab are used to lower and raise indices, respectively. 
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3 The Calculation 

The action during inflation is 



S 



d A x\/—g 



\r + ^ab^(^0) a (^0) b - V(4>) 



(12) 



where g^ v and R are the metric and curvature scalar in the spacetime, and h a b is the metric 
in the scalar field space. The equation of motion for <fi in the spatially flat homogeneous and 
isotropic background is 

a + 3H<j) a + V' a = 0, (13) 

where a = D<fi a /dt = b Vt>0 a and V a is the covariant derivative in the scalar field space. 
We also have 

(14) 



3H 2 = V + ^ a a and H 



2^ 



We define the slow-roll parameters 



and 



2 H 2 

and make the slow-roll assumptions 

e = 0(0 and 
for some small parameter £. We then have 



H\</>\* 
5 = O(0 , 



H 



2e(e + 6) = O^ 2 ) 



(15) 



(16) 



(17) 



and we will make the standard extra assumption^ 



5 

H 



o(e 



3.1 Equation of motion for the perturbations 

The equation of motion for the inflaton perturbation 5<f) a (k, t) on flat hypersurfaces is 



D 



+ 



k 

a 



+ V 



ab i 



1 D (a 3 ■ ■ 



i 3 dt\H 



dt 2 dt 

where R a bcd is the curvature tensor in the scalar field space, defined by 

^bcd^ = (VcV d - V d Ve) V*. 

Defining (p a = a5<f) a , the conformal time drj = dt/a and x = —krj we get 



a b U0 b , (19) 



(20) 



2,„a 



D*<p 



dx 2 



(f) 



21 



aH 



H 

IT 2 



h ah + R a J 



0C0d ^;ab 1 D 



3 HHdt 



a 3 H- 



a 



b a bN 
(21) 



For a more general slow-roll approximation which doesn't make this extra assumption, see Ref. p 



V?b 
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Using 

the equation of motion is 



. dt k 
x = -k I — = —7 
a ah 



1 + e + 3e 2 + 2e5 + O (£ 



whereQ 



-ab 



,2 § A ,ab A 5 \ a b ID (4> a ^ 
+ i(l + 2e)ir cd b |^-(l + 2 £ )^ + 0(f 



3.2 Quantization and Green's function solution 

The quantization conditions are 

2<5 (3) (x-y) /i ab , 



' a/ \ Dip h . , 



9?7 

otherwise zero. With the Fourier transformation 



(2tt) 3 / 2 

the quantization conditions for the Fourier modes are 



ikx 



bt 



^,^,^-(1,7?) 



A 3) (k-l)/i 



ab 



otherwise zero. 

Now, the solution of the homogeneous part of Eq. 

D 2 Vo a + A 2 



dx 2 



with the vacuum boundary condition at x — > 00 is 

1 



99 a (k, x) 



a a (k) ^o(ar) + a ta (-k) <p*\ 



x) 



where 
and 



(k),a tb (l)l =6 (3 \k-\)h 



ab 



(p (x) 



1 H — I e" 



f22) 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 

(29) 

(30) 
(31) 



The Green's function solution of Eq. (^) with these boundary conditions is 

^ POO (£y^ 

<p*(k, x) = <^ a (k, x) + — — C a b (fc, w) y? b (k, u) [cp* (u) Vo {x) - cp* (x) <p (u)) . (32) 

2 Jx u 



4^ab rec j uces ^ e ab Q £ p^ e f g or( [ er £ ; an( j to g/3 of Ref. 0] in the case of a single component inflaton. 



3.3 Slow-roll expansion 

To implement the standardf] slow-roll expansion, we expand ( ah in a power series in In x, 



n=0 



T) ! 



(33) 



and, in addition to Eqs. (|T6| ) and (|18D, assume 

Cn ah = o(n- (34) 

Therefore, to obtain the solution up to second-order corrections, it is sufficient to consider 



c ab = cr+( 2 aD ini 



ab i /■ ab 



where 



ab 



-ab 



-ab 



2 = 1 



aH=k 



(35) 



(36) 



1 - If + -5 } eh ab 4 f I - -J 



ka J,b 



+ 



1 D 



ka Ab s 



1 

3 



3 ) H H 3 Hdt \H H 

c , ji \ r ;ab 



l + 2 e )i? a cd b ^-(l + 2 e) ^ + 0(e 3 ) 



(37) 



and 



ab 



DC 



ab 



din a; 



DC 



ab 



-2e (e + 6) h 



ab 



D 



aH=k 



(38) 



ia ib 1 



HdtXHH 



1 L> 



Ic icT 



3M 



+ 



D V 



ab v 



na b< 

. cd H H I ' Hdt \ 3H 2 



(39) 



aH=k 



5 For a more general slow-roll expansion, see Ref. jl 
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3.4 Power spectrum 

Now, to calculate the power spectrum, we substitute Eq. ( pop into Eq. ( p2|) and integrate 
iteratively. Most of the calculations have already been done in 0, and, in the limit x —>■ 0, 
the asymptotic form for (p a up to second-order corrections is 



¥? a (k,x) 



1 i 



2k x 



a a (k) - a ta (-k) + Z a h {k, x) a b (k) - Z a h *{k, x) a tb (-k) 



where 



7 ah 



ab 



7T 



OC 



—a — 4 H — - + if a 



/■ a /- cb 
SI cU 



+ 



7T" 



a H — a; h vk a + 



12 



C: 



ab 



3 2 

+ ^(Ci a cCi cb -C 2 ab )(lna;) 2 



ab 



lnx 



1 



The numerical constant a is defined as 

a = 2-ln2-7~ 0.729637, 
where 7 ~ 0.577216 is the Euler-Mascheroni constant. Eq. (|40"D gives 



<^ a (k, x) <^ bt (l,x 



2kx : 



<^(k-l) 



^ab + ^ab + ^ba* + /^ac^bd* 



2A;ar 



// ab : 2«Ci ab + ( 2c? - %a - 4 + ^ ) Ci a c Ci cb + f« 2 + ?a - ^ 1 C: 



ab 



12 / 



Ci ab + ( 2a - ± ) Ci a c Ci cb + 



a /- cb 



ab 



and so, from Eq. ([?]), the power spectrum is 



(40) 



(41) 



(42) 



lnx+(2Ci a c Ci cb -C 2 ab )(lnx) 2 } (43) 



V(k) = N a N b 



h ah + 2aCi aD + I 2a z - ~-a - 4 + — I Ci a c Ci CD + [^ + ^a- — \£ 



2ixax 



ab 



X 



7T 



71 



ab 



12 



- 2 



ab 



2a 



1 



Ci a c Ci cb + 3C2 



ab 



lnx+(2d a c Ci cb -C2 ab ) (hi*) 2 } (44) 



The right hand side is constant to order £ 2 and so we can choose to evaluate it at any 
convenient time around horizon crossing. We will evaluate it at aH = k. 
Now, from Eq. 



x\ aH=k ~ 1 + e + 3e 2 + 2e5\ aH -_ 



(45) 
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Therefore 



V(k) 



j- ) N a N h { (l - 2e - 3e 2 - Ae5) h ab + [2a - (4a + 2)e] Ci 



f- 

V2tt 

+ H-I tt - 4+ y)«* + ( tt2 + f a -S) & 



ab 



ab 



aH=k 



where Ci ab and C 2 ab are given by Eqs. fl3"6| ) and 

To write V(k) in terms of the inflaton potential, we define 



_ V> a V a __ V ah 



V a bc = U ^ 



the unit covectors 



the unit vectors 



and the component notation 



N> h N h 



and e 



v 



h ab e£ and e^ = h ab e 



ab V 
b 5 



Now 



wjq = efjWa and Wy = e v w 8 



V 



and using Eqs. (|H|) and ( H ) we can derive 



H 2 = j [l + ~U - hj 2 + l -UV vv + O (C 

e= l -U (l - 2 -U + ^Vvv) + O (?) , 
5 = l -U - V vv + O (e) , 



D 

Hdi 



-W*Vy - \u^h ah V aV Vb 



(46) 



(47) 

(48) 

(49) 

(50) 
(51) 

(52) 

(53) 
(54) 

(55) 
(56) 



The second term in Eq. ([56]) is needed because, for example, V a vb = V a by — UR a wb and 
UR&vvh is of order £. Thus, while V a by is of order £ 2 , V a vb is of order £. We could drop the 
second term in Eq. ( j5"6f) if we assumed either U = 0(^ 2 ), which is physically very reasonable 
but non-standard, or R a wb = which is possible but not what one would expect in 

general. 
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Expressing Eqs. ([46f) , ( 36|) and (|38|) in terms of the potential gives 



2n - [ -Qt + 1] W 



tl-NN 



+ \2a 2 - -a - 4 



7T 



^CliVaClbjV + \ ® A + - ^ J C2JVJV 



7T 



(57) 



where 



Gab = i(l + |w-2V yy )w/ iab +(l + ^)(Were b / -V ab ) 



+ 2^1+ gW^j UR aV Vb + g^^ Cd (-^aVcb + -RacVb) VdV 



(5f 



and 



C 2ab = W (Vyy - W) /lab - 2W 2 eK + WV ab + U (e^Vw + e^V aV ) - V ahV - ^ cd V ab cV d y 

-U 2 R a vVb + g^ Cd (-^aVcb + -RacVb) VdV + g W 2 ^-RaVVb;V + g ^ Cd -RaVVb;cVdy^ • 



(59) 



Substituting Eqs. (|58|) and (|59f) into Eq. fl57|), we can obtain the power spectrum in terms 
of the inflaton potential and its derivatives explicitly. The result is 



V(k) 



VN, C N'' 
12tt 2 



1 



17 5tt 2 



1 + ( a - - ) U + I --a 2 + a L + 

! 6/ V 2 6 24 



U 1 



+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



( 2 § 7tt 2 
2a + 2a 2 - -a - 9 + — 
\ •) o 

-2a+ (-a 2 + + 5 

2 



lAh 2 NV 



4 7tt 2 ^ 



3 

5 2tt 2S 



12 



U 



V 



NN 



u 



UR 



NVVN 



a + —a h 

3 V 39 

2 4 13 7T 2 \ , 

-2a 2 + + 8 - J Uh NV V NV 
3 o / 



2 Z 7T" 

—a a H 

3 12 



1 

— i 

3 



Vatatv + -/i ab ViVAraVby 



77 n + — — \ U 2 ( RnVVN;V + -^ ab -RAfVVAf;aVbV 



1 

— / 

3 



2a A 



7T 



« - 4 + — U ab ( V aN - ^UR NVVa ) I V hN - ^UR N vvb 



1 



2 „ 8 7T 2 



+ Uh ab R NVaN V hV + O (e 



(60) 



Finally, the spectral index, n — 1 = dlnV/dlnk, and the running of the spectral index 
can be calculated from Eq. ( f44|) in a similar way. The results are 

11 
~6 



and 



n-l = -U+[2a- — )U 2 + (-2+U)Uh 2 NV 



+ 



2+ [-2a+-)U 



Vnn + 



h -a + - \ U 

3 V3 9 



UR 



NVVN 



+ ( -2a + ^ ) UVvv + ( 4a - ^ ) Uh NV V NV 



( 1 - - 

+ ^2a + (Vjvjvv + g/i^VjvjvaVby 

(2 2 \ 3 / 1 
— -a — -) U 2 [Rnvvn;V + -R/vwiV;aVbv 

+ (-4a + ^ /i ab (v aN - -UR NVVe ) (v hN - ^UR NVVh 
+ (--a - ^) Uh ah R NVaN V hV + O (e 3 ) (61) 



dn 4 

-2W 2 + 2WVivjv - -U 2 R NVVN + 2WV yy - 4Uh NV V? 



d\nk 3 



— 4 (uh 2 NV — Vnn + -^URnvvn 

— 2 (Vnnv + -^ ab VATjv a Vbv 

2 3 / lb 

+ -W 2 ( Rnvvn-,v + -Rjvyyjv ; aVbv 
+ 4/i ab f V aW - WJ fVbiv - \UR 



3 ""7 V 3 



NVVb 



+ ^Uh ah R NVa NV hV + O (e) ■ (62) 
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